Preprint typeset in .IHEP stylo - HYPER VERSION 



New nonuniform black string solutions 



Burkhard Kleihaus^* Jutta Kunz^^ and Eugen Radu 



'4 



^ Institut fiir Physik, Universitdt Oldenburg, Postfach 2503 D-26111 Oldenburg, Germany 
"^Department of Mathematical Physics, National University of Ireland, Maynooth, Ireland 



Abstract: We present nonuniform vacuum black strings in five and six spacetime dimen- 
sions. Tlie conserved cliarges and the action of these solutions are computed by employing a 
quasilocal formalism. We find qualitative agreement of the physical properties of nonuniform 
black strings in five and six dimensions. Our results offer further evidence that the black hole 
and the black string branches merge at a topology changing transition. We generate black 
string solutions of the Einstein-Maxwell-dilaton theory by using a Harrison transformation. 
We argue that the basic features of these solutions can be derived from those of the vacuum 
black string configurations. 



Keywords: black strings, numerical solutions, Harrison transformation 



*E-mail: kleihausOtheorie .physik.uni-oldenburg.de 
^E-mail: kunzOtheorie . physik . uni-oldenburg . de 
■''E-mail: raduOthphys . nuim . ie 



Contents 



|l|. Introduction 



I 



General ansatz and properties of the solutions 



2.1 The equations and boundary conditions 



2.2 Properties of the solutions 



2.3 A counterterm approach 

2.4 Remarks on nonuniform bubble-like solutions 



Numerical nonuniform black string solutions 



3.1 Numerical procedure 



3.2 Black string properties 



3.3 Black strings and black holes 



^. Black string solutions in Einstein-Maxwell-dilaton theory 

0| Asymptotically M.^^^ x EMD black string solutions 
|4.2| Black strings in a background magnetic field 

^. Conclusions 

Appendix 







16 

20 

21 
23 

25 



1. Introduction 

In recent years it has been realized that, even at the classical level, gravity exhibits a much 
richer structure in higher dimensions than in four dimensions. 

Black string solutions, present for D > 5 spacetime dimensions, are of particular interest, 
since they exhibit new features that have no analogue in the black hole case. Such configu- 
rations are important if one supposes the existence of extradimensions in the universe, which 
are likely to be compact and described by a Kaluza-Klein (KK) theory. 

The simplest vacuum static solution of this type is found by assuming translational 
symmetry along the extracoordinate direction, and corresponds to a uniform black string 
with horizon topology S^^^ x S^. Although this solution exists for all values of the mass, it 
is unstable below a critical value as shown by Gregory and Laflamme Q. This instability was 
interpreted to mean that a light uniform string decays to a black hole since that has higher 
entropy. However, Horowitz and Maeda |^ argued that the horizon could not pinch off, so 
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the end state of the instabihty could not be a cohection of separate black holes. Instead, they 
conjectured that the solution would settle down to a non-translationally invariant solution 
with the same horizon topology as the original configuration. 

This prompted a search for this missing link, and a branch of static nonuniform black 
string solutions was subsequently found by perturbing the uniform black string by the thresh- 
old unstable mode. The D = 5 approximate solutions were obtained in j^, perturbed D = 6 
were presented in Q, and higher dimensional generalizations were discussed in 

Despite some attempts, no analytic solutions are available for nonuniform black string 
solutions and one has to employ numerical techniques. At present, the nonuniform branch is 
numerically known only for D = 6 Q (see also the post-analysis in Q). 

Apart from the black string solutions, KK theory possesses also a branch of black hole 
solutions with an event horizon of topology The numerical results presented in 

(following a conjecture put forward in Q) suggest that, at least for D = 6, the black hole and 
the nonuniform string branches merge at a topology changing transition. Still, a number of 
aspects remain to be clarified, and the literature on nonuniform black string and black hole 
solutions in KK theory is continuously growing (see |lO| , pT| ] for recent reviews). 

The main purpose of this paper is to numerically construct and study the nonuniform 
black string branch in D = 5. This dimension is of particular interest since one may join the 



black string results with those of the D = 5 black hole branch discussed in (see also |12|). 

We begin with a presentation of the general ansatz and the relevant quantities for an 
arbitrary spacetime dimension D. In this context, we propose to compute the mass, tension 
and action of the nonuniform black string solutions by using a quasilocal formalism. In Section 
3 we present our numerical results. We demonstrate that for D = 5 a branch of nonuniform 
black string solutions exists, at least within the scope of our numerical approximation. The 
numerical methods used here are rather different from the methods employed to obtain the 
D = 6 solutions Q. We construct nonuniform black string solutions also in D = 6 dimensions, 
extending the known set of solutions |Q to larger deformation of the event horizon. In 
Section 4 solutions of the Einstein-Maxwell-dilaton (EMD) equations are generated from the 
vacuum configurations, by using a Harrison transformation originally derived in ||l3|. The 
basic properties of two different types of solutions, corresponding to charged black strings 
which asymptote to x and black strings in a background magnetic field are discussed 

there. We give our conclusions and remarks in the final section. The Appendix contains a 
brief discussion of some numerical aspects of our solutions. 

2. General ansatz and properties of the solutions 

2.1 The equations and boundary conditions 

We consider the Einstein action 

I=--J-—[ d^x^/^R--^l d^-^xV^K, (2.1) 



levrG Jm 8ttG jqm 
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in a D— dimensional spacetime. The last term in is the Hawking-Gibbons surface term 
[^], which is required in order to have a well-defined variational principle. K is the trace of 
the extrinsic curvature for the boundary ^^A and h is the induced metric of the boundary. 

We consider black string solutions approaching asymptotically the D — 1 dimensional 
Minkowski-space times a circle Ai^^^ x . We denote the compact direction as z = x^~^ 
and the directions of R ^ ^, while = t. The direction z is periodic with 

period L. We also define the radial coordinate r by = (x^)^ + • • • + 

The nonuniform black string solutions presented in this paper are found within the metric 



ansatz ^ 



where 



/ = l-(^)^-^ 

r 

The Einstein equations = 0, G^. + G^ = and Gg = (where 9 denotes an angle of 
the D — 3 dimensional sphere) then yield for the functions A, B, C the set of equations 

.// A ,,2 J AC A' f'C f \ f" 3fA' 

^'' + -^+^'^ + y + (i^-3) + — + - + + +Z_ + _L_=o, 



(2.3) 



2B-2C 



/' Ac + f'C' A' ^ 

+ —, + 7^ + — + A'C' = 0, 

rf f r 

where a prime denotes d/dr^ and a dot d/dz. 

The remaining Einstein equations = 0, GJ^ — = yield two constraints. Following 
[|, we note that setting G^ = G^ = G^^ = G^ + Gf = in V ^G^''^ = and V^G^^ = 0, we 
obtain 

d. {^gGl) + (y7^i(G; - G^) = 0, (2.4) 

V75. (^g:) - 9, (y7v^l(G; - gd) = 0, 



^An ansatz involving only two undetermined functions was suggested in for a special coordinate choice. 
However, we could not obtain numerical solutions of the Einstein equations within that reduced ansatz. 
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and, defining r via d/df = \/fd/dr, then yields the Cauchy-Riemann relations 




Thus the weighted constraints satisfy Laplace equations, and the constraints are fulfilled, 
when one of them is satisfied on the boundary and the other at a single point Q. 

The event horizon resides at a surface of constant radial coordinate r = tq and is char- 
acterized by the condition /(rg) = 0. Introducing the coordinate f, where r = ^/tq + f^, the 
horizon resides at f = 0. 

Utilizing the reflection symmetry of the nonuniform black strings w.r.t. z = L/2, the 
solutions are constructed subject to the following set of boundary conditions 

A\,^^ = B\,^^ = C\,^^ = 0, (2.6) 



A\,^o - S|,=o = do, drA\,^o = drC\,^o = 0> (2-7) 

dzAUo,L/2 = d.Bl=o,L/2 = dzCUo,L/2 = 0, (2.8) 

where the constant do is related to the Hawking temperature of the solutions. Regularity 
further requires that the condition dfB\-^Q = holds for the solutions. The boundary 
conditions guarantee, that the constraints are satisfied, since ^/^^G\ = everywhere on the 
boundary, and \^\/—g{G'^ — G^) = on the horizon. 

2.2 Properties of the solutions 

For any static spacetime which is asymptotically Ai^^^ x one can define a mass M and a 
tension T [O], these quantities being encoded in the asymptotics of the metric potentials. As 



discussed in [17|, |18|, the asymptotic form of the relevant metric components of any static 
solution is 

9tt^-^ + ^, 9zz^l + ^- (2.9) 

When computing M, T or the gravitational action, the essential idea is to consider the 
asymptotic values of the gravitational field far away from the black string and to compare 
them with those corresponding to a gravitational field in the absence of the black string. 
Therefore, this prescription provides results that are relative to the choice of a reference 
background, whose obvious choice in our case is A4^^^ x S^. 

The mass and tension of black string solutions as computed in ^] are given by 



-4- 



where JId-s is the area of the unit S^~^ sphere. The corresponding quantities of the uniform 
string solution Mq and Tq are obtained from ( 2.10D for = 0, q = r^'^. One can also define 
a relative tension n (also called the relative binding energy or scalar charge) 

n = (2.11) 

which measures how large the tension is relative to the mass. This dimensionless quantity 
is bounded, < n < D — 3. Uniform string solutions have relative tension no = l/(-C — 3). 
Another useful quantity is the rescaled dimensionless mass 

The Hawking temperature and entropy of the black string solutions are given by 

= gAo-Boj.^^ S = Sq\ re^°+(^-3)Coj^^ (2.13) 
L Jo 

where Tq, Sq are the corresponding quantities of the uniform solution 



So = -mD-sr^-^ (2.14) 



and Ao{z), Bo{z), Co{z) are the values of the metric functions on the event horizon r = rg. 
Together with the mass M and relative tension n, these quantities obey the Smarr formula 

[0] 

TS = ^^^~" m. (2.15) 

Note that the relations ( |2.9D -( 2T^ ) and ( p.l5| ) are also valid for black hole solutions. 

Black string thermodynamics can be discussed by employing the very general connec- 
tion between entropy and geometry established in the Euclidean path integral approach to 
quantum gravity |14|. In this approach, the partition function for the gravitational field is 
defined by a sum over all smooth Euclidean geometries which are periodic with a period /3 
in imaginary time. This integral is computed by using the saddle point approaximation, and 
the energy and entropy of the solutions are evaluated by standard thermodynamic formulae. 

We consider a canonical ensemble with Helmholz free energy (thus at fixed temperature 
and extradimension length) 

F[T,L] = ^ = M -TS (2.16) 

The Euclidean action of the vacuum solutions computed by subtracting the background con- 
tribution is 

I = Y^^^D-3ict - c.) = ^^{M + TL), (2.17) 
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with /3 = l/T. 

The first law of thermodynamics reads 



dM = TdS + TdL. (2.18) 

It follows that 



Combining the Smarr formula ( 2.15 ) and the first law, it follows that, given a curve n{ji) in 



the (n, ;^)-plane, the entire thermodynamics can be obtained |17]. 

We remark also that the Einstein equations (|2.3|) are left invariant by the transformation 
r — > fcr, z kz, — > ro/k, with k an arbitrary positive integer. Therefore, one may generate 
a family of vacuum solutions in this way, termed copies of solutions. The new solutions have 
the same length of the extradimension. Their relevant properties, expressed in terms of the 
corresponding properties of the initial solution, read 

= -j^Zi^ Tk = kT, Sk = -j^^^ "fc = (2-20) 
This transformation, suggested first by Horowitz [^], has been discussed in [^] for D = 6. 
2.3 A counterterm approach 

Similar results for the black strings' mass, tension and action are obtained by using the 



quasilocal tensor of Brown and York [22|, augmented by the counterterms formalism. This 
technique consists in adding suitable counterterms Id to the action of the theory. These 
counterterms are built up with curvature invariants of the induced metric on the boundary 
dM (which is sent to infinity after the integration) and thus obviously they do not alter the 
bulk equations of motion. By choosing appropriate counterterms which cancel the divergen- 
cies, one can then obtain well-defined expressions for the action and the energy momentum 
of the spacetime. Unlike the background substraction, this procedure is satisfying since it is 
intrinsic to the spacetime of interest and it is unambiguous once the counterterm is specified. 
While there is a general algorithm to generate the counterterms for asymptotically (anti-)de 



Sitter spacetimes, the asymptotically fiat case is less-explored (see however |23] and the more 



general approach in recent work |24|). 

It is also important to note that the counterterm method gives results that are equivalent 
to what one obtains using the background subtraction method. However, we employ it because 
it appears to be a more general technique than background subtraction, and it is interesting 
to explore the range of problems to which it applies. 

Therefore, we add the following counterterm part to the action principle (p.lj) 
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where TZ is the Ricci scalar of the boundary geometry. 

Varying the total action with respect to the boundary metric hij, we compute the bound- 
ary stress-tensor 

Tij = ^T^n = ^ (^ij - h^^K - - nh,,) - h,,u^ + , (2.22) 



-h 6h^^ SvrG 



where Kij is the extrinsic curvature of the boundary and ^' = y (,0^4)7^ • 

Provided the boundary geometry has an isometry generated by a Killing vector a 
conserved charge 

Qi = ^ d'^-^S' ^Tij (2.23) 

can be associated with a closed surface S. Physically, this means that a collection of observers 
on the boundary with the induced metric hij measure the same value of Q^. 

The mass and tension are the charges associated to d/dt and d/dz, respectively (note 
that d/dz is a Killing symmetry of the boundary metric). The relevant components of the 
boundary stress tensor are 



The mass and tension computed from ( |2.23| ) agree with the expressions ( p. 10 ) 



2 



One should remark that the counterterm choice is not unique, other choices being possible 
as well (see |25] for a related discussion). Our choice of using ( 2.21|) was motivated by 



the fact that the general expression for the boundary stress-tensor is very simple. For an 
asymptotically Ai^^^xS^ spacetime, we find that Id (2.21) or other possible choices proposed 



in ||2^ regularizes also the black strings Euclidean action, yielding similar results to those 
obtained within the background subtraction approach. 

2.4 Remarks on nonuniform bubble-like solutions 

There is a also a simple way to generate time dependent solutions with a nontrivial extradi- 
mension dependence. A Lorentzian solution of the vacuum Einstein equations is found by 
using the following analytic continuation in the general configuration ( |2.2| ) (with dQ^_^ = 

de^ + sin^ ednl_^) 

t = ix, 9-7^/2 = IT. (2.25) 

The new solution reads 

ds^ = e2^(^'^)/(r)dx' + e^^^'-'^) f ^ + dzA + e^^^^-'^^r^-dT^ + cosh^ rdQl.,), (2.26) 



^Note that in computing T, one should consider the integration over 5*^ — /At, with At — f dt. However, 
a similar problem appears also in the Hamiltonian approach jl9|]. 



-7- 



with A{r,z), B(r,z), C{r,z) and /(r) the functions of the black string configuration. 

This technique of double analytic continuation was originally developed for the study of 
the stability of the KK vacuum [^], ||2^, and has been considered in the last years by many 
authors in AdS/CFT context (see e.g. 

However, the solution one finds starting with a nonuniform vacuum black string has some 
new features, as compared to the Schwarzschild-Tangherlini seed solution. In that case, the 
resulting spacetime describes a contracting and then expanding "bubble of nothing". The 
new configuration ( |2.26| ) has now two compact extradimensions z and x- The radial variable 
r is restricted to the range r > tq. r = tq is not the boundary of spacetime, but it is the 
S^~^ X of minimal area. The curves at r = tq with constant z (and constant points on 
S^~^) are geodesies. Also, regularity at r = tq requires x to be periodic with period (3 = 1/T, 
with T given by ( p. 13 ). One can see that the geometry traced out by the r = vq surface is 



the (D — 3)-dimensional de Sitter spacetime with a 2— dependent conformal factor, times the 
extradimension. Thus ( 2.26| ) would rather describe a nonuniform "vortex of nothing". 



One can use the same techniques as in Section 4 to find the corresponding D— dimensional 
solutions in Einstein-Maxwell-dilaton theory. 

KK bubble solutions in d = 5, 6 dimensions have been also considered in Ref. [p9|| . A 
number of exact solutions have been presented there, describing sequences of KK bubbles 
and black holes, placed alternately so that the black holes are held apart by the bubbles. 
However, the configurations in ||2^ differ from ( 2.26| ), since they are static and asymptotically 
approach Minkowski spacetime times a circle. 

3. Numerical nonuniform black string solutions 
3.1 Numerical procedure 

Our main concern here is the numerical construction of nonuniform black string solutions in 
D = 5 dimensions. We have also constructed nonuniform black string solutions in I? = 6 
dimensions, reproducing and extending the known set of solutions, obtained previously by 
different methods, which could not be successfully applied in D = 5 dimensions ^, 0, 

To obtain nonuniform black string solutions, we solve the set of three coupled non- 
linear elliptic partial differential equations numerically |30|, subject to the above boundary 
conditions. We employ dimensionless coordinates f and z, 

r = f/{l + r), z = z/L, (3.1) 

where the compactified radial coordinate f maps spatial infinity to the finite value f = 1, and 
L is the asymptotic length of the compact direction. The numerical calculations are based 



on the Newton-Raphson method and are performed with help of the program FIDISOL |3C], 
which provides also an error estimate for each unknown function. 

The equations are discretized on a non-equidistant grid in f and z. Typical grids used 
have sizes 65 x 50, covering the integration region < f < 1 and < z < 1/2. (See |^ 
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and ||3l| for further details and examples for the numerical procedure.) For the nonuniform 
strings the estimated relative errors range from approximately ~ 0.001% for small geometric 
deformation to ~ 1% for large deformation. Further discussion of the numerical accuracy is 
deferred to Appendix A. 

The horizon coordinate ro and the asymptotic length L of the compact direction enter 
the equations of motion as parameters. The results presented are mainly obtained with the 
parameter choice 

ro = l, L = L-' = |"^" ° = ^ (3,2) 
" ' I 4.9516 D = 6 ' ^ ' 

where L*^"* represents the value, where the instability of the uniform string occurs. The 
branch of nonuniform strings is then obtained by starting at the critical point of the uniform 



string branch and varying the boundary parameter do, which enters the Eq. (2/7), relating 
the values of the functions A and B at the horizon. 

3.2 Black string properties 

Let us first consider the metric functions A, B and C for nonuniform string solutions as 
functions of the radial coordinate r and of the coordinate z of the compact direction. Keeping 
the asymptotic length L of the compact direction and the horizon coordinate rg fixed, the 
solutions change smoothly with boundary parameter d^. A measure of the deformation of the 
solutions is given by the nonuniformity parameter A 

X = \( ^ - 1^ , (3.3) 



2 V7^n____ 

where 7?.max and T^min represent the maximum radius of a (D — 3)-sphere on the horizon and 
the minimum radius, being the radius of the 'waist'. Thus for uniform black strings A = 0, 
while the conjectured horizon topology changing transition should be approached for A ^ co 
^ . As do first increases and then decreases again, the nonuniformity parameter A increases 
monotonically. 

In Figure 1 we exhibit the metric functions for D = 5 nonuniform string solutions for 
several values of the nonuniformity parameter, A = 1, 2, 5, 9. The functions exhibit extrema 
on the symmetry axis z = at the horizon. As A increases, the extrema increase in height 
and become increasingly sharp (implying a deteriorating numerical accuracy for large values 
of A). 

To obtain a more quantitative picture of the metric functions, we exhibit —gtt/f = e^"^, 
gzz = e^^, and gee/f'^ = e^*^ in Figure 2 for several fixed values of z. 

The spatial embedding of the horizon into 3-dimensional space is shown in Figure 3 
for the D = 5 nonuniform black string solutions. In these embeddings the proper radius 
of the horizon is plotted against the proper length along the compact direction, yielding a 
geometrical view of the nonuniformity of the solutions. 

The deformation of the horizon of the D = 5 nonuniform black string solutions is further 
explored in Figure 4. The maximum radius of the 2-sphere on the horizon T^max and the 
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X=1.0 



X=2.0 



A(r,z) A(r,z) 




Figure 1. The metric ftmctions A, B and C of the D = 5 nonuniform string solutions are shown as 
fmictions of the radial coordinate r, with the horizon located at rg = 1, and the coordinate z of the 
compact direction with asymptotic length L = L'^'''* = 7.1713, for several values of the nonuniformity 
parameter, A = 1 (first column), A = 2 (second column), A = 5 (third column), A = 9 (forth column). 
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X=5.0 



X=9.0 



A(r,z) A(r,z) 




Figure 1. Figure 1 continued. 
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Figure 2. The metric functions , , and e^*^ of the D = 5 nonuniform string solutions are 
shown as functions of the compactified radial coordinate x — f, for several fixed values of the coor- 
dinate z of the compact direction {z/L = 0, 0.1, 0.2, 0.3, 0.4, 0.5), as well as of the nonuniformity 
parameter A (A = 1 (first column), A = 2 (second column), A = 5 (third column), A = 9 (fourth 
column)). 
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Figure 2. Figure 2 continued. 
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Figure 3. The spatial embedding of the horizon of _D = 5 nonuniform black string solutions with 
horizon coordinate ro = 1 and asymptotic length of the compact direction L = L'^"* = 7.1713, is 
shown for several values of the nonuniformity parameter, A = 0.5, 1, 2, 3, 5, 9. 
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minimum radius T^min representing the 'waist' are presented together with the proper length 
Lh of the horizon along the compact direction as functions of the nonuniformity parameter 
A, ranging from < A < 9. 

With increasing A, T^max first increases, reaches a maximum around A ~ 4 and then 
decreases slightly again towards still larger values of the nonuniformity parameter; in contrast 
Ln/L increases monotonically and T^min decreases monotonically (in the range considered). 
We expect that T^max and Ln/L approach finite values in the limit A — > oo, whereas T^min 
should reach zero in this limit, when extrapolated (approximately linearly) in the figure. For 
comparison, we also show in the figure the corresponding geometric quantities of the D = 6 
nonuniform black string solutions. 




Figure 4. The maximum radius of the (D — 3)-sphere on the horizon T^-max, the minimum radius 7?.niin, 
and the proper length Ln of the horizon along the compact direction divided by the asymptotic length 
L — L"^^ are shown for D ^ 5 and D = 6 nonuniform black string solutions as functions of 1 / ( 1 + A) . 



We exhibit in Figure 5 the mass M, the relative tension n, the temperature T and the 
entropy S of the D = 5 and D = 6 nonuniform string solutions, in units of the corresponding 
uniform string solution, versus the nonuniformity parameter A. Interestingly, the mass and 
the entropy assume a maximal value in the vicinity of A ~ 4, while the tension and the 
temperature assume a minimal value, both in 5 and in 6 dimensions. Since the extrema 
appear only around A ~ 4, which is the maximal value of A obtained in previous calculations 
[^, they were not recognized there. Extrapolating these quantities to A — > cxd yields for the 
tension the critical value n*, where n*/no ~ 0.8 for D = 5 and n^/riQ ^ 0.6 in six spacetime 
dimensions. 

The mass and tension exhibited in Figure 5 are obtained from the 1st law of thermo- 
dynamics together with the Smarr relation ( |2.I5| ). A discussion of the mass and the string 
tension as obtained from the asymptotic fall-off of the metric functions is given in Appendix 
A. 
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Figure 5. The Af, the relative tension n, the temperature T and the entropy S of the D = 5 
(a) and D = 6 (b) nonuniform string solutions are shown in units of the corresponding uniform string 
solution (denoted by Mq, no, Tq, and So) as functions of 1/(1 + A). 



3.3 Black strings and black holes 

In I? = 6 dimensions, evidence was provided that the nonuniform string branch and the black 
hole branch merge at a topology changing solution |^]. We would now like to reconsider this 
evidence in the light of the continuation of the D = Q nonuniform string branch to larger 
deformations, and further address the question, whether there is analogous evidence in D = 5 
dimensions. 

We therefore exhibit in Figure 6 the mass M versus the relative string tension n, for 
the nonuniform string branch as well as for the black hole branch in 5 and 6 dimensions. 
The black hole data are taken from [^. First of all we note qualitative agreement of the 
shape and the relative position of the nonuniform string branch and the black hole branch 
in 5 dimensions with the shape and relative position of the corresponding branches in 6 
dimensions. But compared to the data and discussion given in Q, we here observe a new 
feature: the backbending of the nonuniform string branch at a critical (minimal) value of the 
relative string tension n^. Although the onset of this backbending can already be anticipated 
in the D = 6 data of But the backbending of the nonuniform string branch at rif, is still 
in accordance with the conjecture of a topology changing transition, occurring at > nh, 
both in 5 and 6 dimensions. 

We attribute the presence of the gap between the black hole branch and the nonuniform 
string branch mainly to insufficient numerical data of the black hole branch close to the 
anticipated transition point. At such a transition point the nonuniform string parameter 
i^min must approach zero (see Figure 4), and likewise the black hole parameter Laxis must 
approach zero, where Laxis measures the proper length along the exposed symmetry axis 
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Figure 6. The mass M of the D = 5 (a) and D = 6 (b) nonuniform string and black hole branches is 
shown versus the relative string tension n. M and n are normalized by the values of the correspond- 
ing uniform string solutions. Here and in Figures 7-11, the data for the black hole branches is from B . 



We note, that close to the anticipated transition point, i?min has decreased considerably 
farther (at the last numerically obtained point of the nonuniform string branch), than Laxis 
(at the last numerically obtained point of the black hole branch). 

Since Laxis has decreased farther for D = Q black holes than for D = 5 black holes, the 
gap between the branches is smaller in 6 dimensions than in 5 dimensions. We exhibit -RmiiD 
-Rmax) and Laxis Eis Well as the black hole equatorial radius i?eq in Figure 7 for D = 5 and 
D = 6 solutions. i?min and Rmax both exhibit the backbending feature present for nonuniform 
string solutions at large deformations. The figure is consistent with the vanishing of Laxis and 
-^min at the same critical value of n. There R^q and iJmax should also merge. The transition 
might then occur in the vicinity of n^,/nQ ~ 0.8 for D = 5 and n^^/riQ ~ 0.6 for D = Q (as 
opposed to nb/no ~ 0.55 for D = 6, which was earlier assumed to be the transition point [^, 
but which is now realized to be the point where the backbending occurs). Extrapolating the 
black hole branch towards this critical value, the i?eq curve appears to smoothly reach the 
endpoint of the (backbending) upper part of the i?max curve of the nonuniform string branch. 
The black hole data are again taken from Q. 

Addressing the thermodynamic properties of the solutions, we exhibit in Figure 8 and 
Figure 9 the temperature and the entropy oi D = b and D = 6 nonuniform strings and black 
holes. Extrapolating the black hole branch towards the critical value n*, where the transition 
might occur, the black hole curves for temperature and entropy also appear to smoothly 
reach the endpoints of the corresponding (backbending) upper parts of the nonuniform string 
branch. This also holds for the mass, of course. Again, the black hole data are from 
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Figure 7. i?min, .Rmax, ^axis and i?oq of the D = 5 (a) and D = 6 (b) nonuniform string and 
black hole branches are shown versus the relative string tension n. 



For very small masses localized black holes are entropically favoured, and for very large 
masses only uniform strings exist When the entropy is plotted versus the mass for the 
black hole branch and the uniform and nonuniform string branches, one observes, that the 
uniform strings become entropically favoured at a certain value of the mass, lying above the 
critical string mass and below the maximum black hole mass |^. This is illustrated in Figure 
10 for solutions in five and six dimensions. 



D=5 D=6 




0.2 0.4 0.6 O.B 1 0.2 0.4 0.6 0.8 1 
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Figure 8. The temperature T of the D = 5 (a) and D — 6 (h) nonuniform string and black hole 
branches are shown versus the relative string tension n (in units of the uniform string quantities To 
and no). 
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Figure 9. The the entropy S of the D = 5 (a) and L* = 6 (b) nonuniform string and blaclc hole 
branches are shown versus the relative string tension n (in units of the uniform string quantities 
and no). 

Concluding, we observe qualitative agreement of all the physical properties of the solutions in 
5 and in 6 dimensions. This strongly suggests, that the same phenomenon is present in both 
cases. In particular, all data are consistent with the conjecture that the black hole branch and 
the nonuniform string branch merge in a topology changing transition. Our new nonuniform 
string solutions give further credence to this scenario, but they still cannot confirm it. 




Figure 10. The product of entropy and mass SM'^/^ (a) and SM'^^^ (b) for the £» = 5 (a) 
and D = 6 (h) nonuniform string branch, the uniform string branch and the black hole branch is 

shown as a function of the mass (in units of the corresponding critical string quantities). The hori- 
zontal lines represent the curves for the corresponding Schwarzschild-Tangherlini solutions. 
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4. Black string solutions in Einstein-Maxwell-dilaton theory 



We now consider the action describing a gravitating Maxwell field coupled with a dilaton in 
a L)— dimensional spacetime 



(4.1) 



where F = dA. The free parameter a governs the strength of the coupling of the dilaton to 
the Maxwell field. 

The corresponding Einstein-Maxwell-dilaton (EMD) field equations are 

= -^e-^-^F', (4.2) 



with stress-energy tensor 



t},'!) =d^cl)d,^-^g^,\d^\^ (4.3) 

r(r) = e-^'^'^iF.pF^y^ - \g,uF'). 

Finding nonuniform black string and black hole solutions of these equations for a compact 
extradimension constitutes a formidable technical task. However, it appears possible to use 
the 'hidden' symmetris of the model (|4.1| ) arising in the dimensionally reduced theory in order 
to generate nontrivial EMD solutions without actually solving the full set of equations. 

This technique, known for the = 4 Einstein-Maxwell theory since long, was generalized 
in the last years to higher dimensions and various other matter fields. In this paper, we'll 



follow the approach in |13], supposing the existence of one Killing vector d/dy, and write the 
L>— dimensional line element in the form ^ 

ds^ = gyy{x)dy'^ + hij{x)dx^dx\ (4-4) 

performing the KK reduction with respect to the y— direction. As proven in [^], the reduced 
action corresponds to a non-linear cr-model, whose target space possesses a rich geometric 
structure. These symmetries imply the existence of a Harrison transformation nontrivially 
acting on the spacetime variables and matter fields. As a result, one may generate nontrivial 



^Note that the Ref. discusses a more general case, with an action principle containing a (d + 1) 

differential form. Therefore, the general metric ansatz (4.4) will contain a set of d-coordinates y'. 
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solutions of the D— dimensional EMD equations starting with known vacuum configurations. 
A detailed description of this procedure and the explicit form of the Harrison transformation 
is given in Ref. [|l3| (see also the results in |^2|). This approach is valid for any value of 
the dilaton coupling constant (in particular also for Einstein-Maxwell theory, i.e., for a = 
0, = 0), and appears to be different from other results in literature. Ref. |3^], for example, 
"charges up" the neutral KK solutions by uplifting them to eleven dimensional M-theory 
and employing boost and U-duality transformations, which fixes a particular value of the 
dilaton coupling constant a. The resulting solutions of type IIA/B string theory describe 
non-extremal p-branes on a circle. 

Here we present only the resulting solutions, which have rather different properties, de- 
pending whether d/dy is a timelike or spacelike Killing vector. Although the same generation 
techniques apply also for black hole solutions, we'll restrict to the black string case. 

4.1 Asymptotically M^~^ x EMD black string solutions 

We start with a vacuum black string solution, written in the following form 

ds'^ = -V{x)dt'^ + hij{x)dx'dxK (4.5) 

where d/dt is a timelike Killing vector. 

The Harrison transformation in this case generates a one parameter family of black string 
solution in EMD theory, with line element 

ds"^ = -y(cosh2 p - sinh^ /5F)-2"(«-3)^^2 ^ ^^^^^2 ^ _ ^^^^^2 pvf'^hijdx'dx\ (4.6) 
and matter fields 

tanh/3 e"*oy 



(p = (po- 2a(L» - 2)a log(cosh2 p - sinh^ f3V), 

where (f)Q are arbitrary real constants and 

a = {2a'^{D-2) + D-3)-^. (4.8) 

For Einstein-Maxwell theory, we find a = 1/{D — 3), the corresponding value in the Kaluza- 
Klein case being a = 1/{2{D — 2)). 

Both uniform and nonuniform solutions of the EMD theory can be generated in this way. 
For example, the uniform string solution constructed within the metric ansatz (p. 2]) reads 



ds' = ^-^ ,^^dt' + {l + C-^f-Umh'(3f»{^ + dz' + r'dnl_,), 

A, = V2iD^ tanh/?e-^o/(.) 
A* V V ; l + (2i)^-4sinh2/3 ^ ^ 

(p = cj)o- 2a(D - 2)alog(l + (— )^-Sinh2 /3), 

r 
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For all a, the surface r = is an event horizon, while r = is a curvature singularity. The 
extremal limit is found by taking /3 — > oo together with a rescaling of rg and has the form 



\ = ^^^^j^^^fe^^M* > <l> = ^o- 2a{D - 2)a log(l + {^)^-% (4.10) 



c being a real constant. Solutions describing several extremal black strings do also exist |2C]. 

Returning to the general nonuniform string case, we observe that its relevant properties 
can be derived from the corresponding Z)— dimensional vacuum solution. The first thing to 
note is that the causal structure of the region r > tq is similar to the vacuum solutions; in 
particular one finds the same location of the event horizon. For the metric ansatz (|2.2| ), the 
spacetime still approaches the x background as r ^ oo, while the matter fields 

behave asymptotically as 

+ ^^^0 + ^,, (4.11) 

where Qe and Qd correspond, in a suitable normalization, to the electric and the dilaton 
charges, respectively, $ being the electrostatic potential difference between the event horizon 
and infinity. 



$ = ^/2{D - 2)a e'^'^o tanh/3, 
Qe = -^j2{D-2)a e^'^o sinh/?cosh/3 (4.12) 
Qd = -2aa{D - 2)smli^(3ct. 

The mass M, the string tension T and the relative string tension n of the EMD solutions are 

M = M(l + 2{D-3- n)a sinh^ /?), 

T = T, (4.13) 

n 

n = 7^ — . 

l + 2(L'-3-n)asinh2/? 

The electric charge and the dilaton charge can also be expressed via (note also that the dilaton 
charge is not an independent quantity) 

Qe = e'^<^Osinh2/3, (4.14) 

(^-4) \j2{D-2) ^' ^ ' 

^'^ = -;^3^(^-3-n)Qsinh2/3. 

The relation between Hawking temperature T and the entropy S of the EMD solutions and 
the corresponding quantities T and S of the vacuum seed solution is 

f = r(cOsh /3)-2a(D-2) ^ g ^ g^^^^^ ^)2a(D-2) ^ (4_^5) 
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thus the product TS remains invariant under the Harrison transformation. 

The Smarr relation ( 2.15| ), derived in |17] for the vacuum case, admits a straightforward 



generahzation to EMD theory, 

^^l^M^TS-^-^^MR^^Q^, (4.16) 

D-2 D-2 ^ ^ ^ 

where Qe = ^D-?,LQe- Therefore the thermodynamics of the EMD solutions can be derived 
from the vacuum solutions. When the parameter (5 is large one has a near extremal charged 
black string. However, a discussion of the extremal limit seems to require knowledge of the 
region r < tq of the seed metric. 

We conclude that every vacuum solution is associated with a family of charged solu- 
tions, which depends on the parameter (3. In particular, the branch of non-uniform solutions 
emerging from the uniform black string at the threshold unstable mode thus must persist for 
strings with non-zero electric charge. The fact that the 'phase diagram' of static solutions is 
qualitatively unchanged as the charge varies strongly suggests that there is still an instability 



for charged black strings |35|. 



Also, a discussion of the thermodynamical properties of these solutions appears possi- 



ble. This is interesting in connection with the Gubser-Mitra conjecture [34|, that correlates 
the dynamical and thermodynamical stability for systems with translational symmetry and 
infinite extent. This conjecture has been discussed by several authors in the last years (see 



e.g. |35]-|Q). Ref. [^] uses the boost/duality transformation to map the phases of KK solu- 
tions onto phases of non- and near-extremal Dp-branes with a circle in their transverse space. 
The results there (see also ||36[] ) confirm the validity of the Gubser-Mitra conjecture for non- 



extremal smeared branes. Similar conclusions are found in Ref. [37], where a discussion of 



thermodynamical stability of charged black p— branes within third-order perturbation theory 
is presented, the Gregory-Laflamme critical wavelength being also determined. 

Asymptotically M^^^ x EMD black hole solutions can be generated by applying the 
same approach, starting with vacuum seed solutions written in the form ([4.5|). The resulting 
solutions are still given by ( |4.6| )-( ]477| ), with the corresponding expressions of V, hij. Similar 
to the black string case, their properties are completely determined by the vacuum seed black 
hole solutions, the Smarr relation ( 4.16| ) being also satisfied. 



4.2 Black strings in a background magnetic field 

A rather different picture is found in the case when d/dy in ( [4.4| ) is a spacelike Killing vector. 
Here we start with a vacuum black string solution ( p. 2D written in the form 

fir) 

sin^ Odif"^ + cos^ edn%_^) 

with 9 G [0, 7r/2], except for D = 5 where £ [0, tt]. 
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After applying a Harrison transformation to this configuration with respect to the Kilhng 
vector d/dy = d/dip, one finds the following solution of the EMD equations 



ds^ = A^'-{r,z,9)( -e^^^'^'^f{r)dt^ + e^''^'^'\4^+dz^) (4.18) 
V fir) 



m 

+e-(-).^(d^^ + cos^ ednU)) + 

^ , ^, 5oe2^("'^)r2sin2 0^ ,^ ^ , 

Mr,z,e) = ^^^—^5,,, (4.19) 

<Pir,z,e) = -^^2(|^^§y^^logA(.,.,^), (4.20) 

where 

i?o is a free parameter which characterizes the central strength of the magnetic field, and a 



is still given by (4. 



One can easily see that this solution is, in terms of the usual definitions, a black string, 
with an event horizon and trapped surfaces. Again, the causal structure of the solution is not 
affected by the Harrison transformation. 



To clarify it's asymptotics, we note that for r — > oo, the line element ( 4.1^ ) approaches 
the simpler form 

ds^ = A2"(-dt2 + dp'^ + dz'^ + dzj + ...+ dzl_^) + A-2"(^-3)^2^^2^ (422) 

A = 1 + , ° 

2a{D-2y ' 

where we introduced the new coordinates {p, zi, . . . zd-a} satisfying rcosO = {zf + . . . + 
dz|)_4)^/^, rsin0 = p, such that [d{rcos6)]'^ + {rcos9)'^dQjy_^ = dzj + . . . + dz]-f_^. 



Therefore, asymptotically the solution (4.18) approaches the Melvin fluxbrane found in 



[pq ] , which represents a higher dimensional generalization of the four dimensional Melvin solu- 
tion. The Melvin magnetic universe is a regular and static, cylindrically symmetric solution to 
Einstein-Maxwell(-dilaton) theory describing a bundle of magnetic flux lines in gravitational- 
magnetostatic equilibrium [^]. This solution has a number of interesting features, providing 
the closest approximation in general relativity for a uniform magnetic field. There exists a 
fairly extensive literature on the properties of this magnetic universe, of particular interest 
being the black hole solutions in universes which are asymptotically Melvin EO]. Black hole 



solutions in a higher dimensional Melvin universe have been constructed recently in |41| 



Therefore, it is natural to interpret the EMD solution ( [4.18 ) as describing a black string 



in an external magnetic field. Note that, similar to the asymptotically A4^ ^ x case, one 
may generate also KK black hole solutions in a Melvin fluxbrane background. 
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One can compute the mass and tension of the black string solutions by using the back- 
ground subtraction approach. In this case, the natural background is the Melvin solution 
( [4.22| ). It follows that, different from the case of asymptotically A4^^^ x solutions, these 
quantities are still given by ( p.lO ). Moreover, it can be proven that the thermodynamic prop- 
erties of the black string EMD solutions are unaffected by the external magnetic field, i.e. 
one finds the same expressions for the Hawking temperature and entropy as in the Bq = no 
magnetic field case. A similar property has been noticed in [^2| for a D = 4 Schwarzschild 
black hole in a Melvin universe background (see also the D > 4 extensions Therefore, 
this seems to be a generic property of uncharged black hole/black string solutions in a back- 
ground magnetic field extending to infinity. Heuristically, this is due to the fact that, in the 
static case, the mass-point/string source of these configurations does not interact directly 
with the background magnetic field. 



5. Conclusions 



Our first concern has been the numerical construction of Z? = 5 nonuniform strings, represent- 
ing solutions of the vacuum Einstein equations. While their physical properties are similar 
to those of -D = 6 nonuniform strings, their construction is more difficult. We attribute this 
to the slower asymptotic fall-off of the metric functions. 

The branch of nonuniform strings emerges smoothly from the uniform string branch at 
the critical point, where its stability changes Keeping the horizon coordinate tq and 
the asymptotic length L of the compact direction fixed, the solutions depend on a single 
parameter, specified via the boundary conditions. Varying this parameter, the nonuniform 
strings become increasingly deformed, as quantified by the nonuniformity parameter A. For 
the largest value of A reached, the 'waist' of the string has a minimal radius of Rmin ~ 0.1 
{Rmin ^ 0.13) for A = 9 (A = 6) in D = 5 {D = 6) dimensions, indicating that it is shrinking 
towards its asymptotic value of zero when A ^ cxd. 

Previously, in D = 6 dimensions, evidence was provided that the nonuniform string 
branch and the black hole branch merge at a topology changing transition [^. Although we 
see a backbending of the nonuniform string branch in both D = 5 and D = 6 dimensions, not 
observed previously, because the nonuniform string branch had not been continued to suffi- 
ciently high deformation, all our data are consistent with the assumption, that the nonuniform 
string branch and the black hole branch merge at such a topology changing transition. In 
fact, extrapolation of the black hole branch towards this transition point appears to match 
well the (extrapolated) endpoint of the (backbending) part of the nonuniform string branch. 

For the phase diagram this would mean that we would have a region < n < nb with 
one branch of black hole solutions, then a region nj-, < n < n,,, with one branch of black 
hole solutions and two branches of nonuniform string solutions, the ordinary one and the 
backbending one, and finally a region < n < uq with only one branch of nonuniform string 
solutions. (We here do not address the bubble-black hole sequences present for n > no). 
Thus the topology changing transition would be associated with n^, and ny, < n < 
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would represent a middle region where three phases would coexist, one black hole and two 
nonuniform strings. This anticipated phase diagram is exhibited in Figure 11. 

This is strongly reminiscent of the phase structure of the rotating black ring-rotating 
black hole system in D = 5 [^3|. The (asymptotically flat) rotating black holes have 
horizon topology, and the (asymptotically flat) rotating black rings have S'^ x horizon 
topology. The rotating black holes exist up to a maximal value of the angular momentum 
(for a given mass) , < J < J* , the rotating black rings are present only above a minimal value 
of the angular momentum (for a given mass), Jb < J, and in the middle region Jb < J < J* 
three phases coexist, one black hole and two black rings p3[ |. 

Further numerical work for nonuniform strings and in particular for black holes in the 
critical region close to n* might confirm this picture further, and it might lead to further 
insight into the structure of the configuration space. The backbending of the nonuniform 
string branch clearly indicates that the configuration space is richer close to the anticipated 
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Figure 11. The mass M of the £) = 6 nonuniform string and black hole branches is shown ver- 
sus the relative string tension n. M and n are normalized by the values of the corresponding uniform 
string solutions. The black hole branch is extrapolated towards the anticipated critical value 7i*. 
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topology changing transition, which currently still invites speculations [0, 44]. 

Our second concern has been the construction of black strings in EMD theory, obtained 
via a Harrison transformation. Different from other results in the literature (see e.g. |33|), 
the construction we proposed in Section 4 is valid for any value of the dilaton coupling 
constant. Apart from asymptotically j\4^~^ x charged black strings, we found also black 
string solutions in an external magnetic field. We found that the properties of these EMD 
configurations can be derived from the corresponding D— dimensional vacuum solution. 

To push forward our understanding of these issues, it would be interesting to find both 
black hole and nonuniform black string branches for D > 6. The dimensions around D = 13 
are of particular interest, since as found in H, for D > 13 the perturbative nonuniform 
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strings are less massive than the uniform solutions. Moreover, their entropy is larger than 
the entropy of uniform strings with the same mass. 
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A. Appendix 

As noted before for black holes, the most error prone part of the numerical procedure is the 
extraction of the mass M and the relative tension n from the asymptotic form of the metric 
via the coefficients ct and Cz @, |6|, ^ . Whereas this extraction is quite accurate for nonuniform 
strings in 6 dimensions, it is quite problematic for nonuniform strings in 5 dimensions. 

To see this we consider the asymptotic expansion for the metric functions in D = 5 
dimensions. 

A- 

and in D = 6 dimensions 

A 

In D = 5 the coefficients must satisfy Coo = ^oo + 2i?ooi a relation essential for the first law 
to hold. 

In 6 dimensions the asymptotic fall-off is sufficiently fast, yielding excellent agreement for 
the mass and tension as obtained from the expansion coefficients with those obtained from 
the first law and the Smarr relation |6|-[p. 

In 5 dimensions the numerical situation is much trickier, because of the presence of the log 
term in the function C. This is aggravated by the fact, that the coefficient Coo is an order of 
magnitude smaller than the other two coefficients. We therefore do not see the log dependence 
in the numerical results, when calculated in the full interval [0 < f < 1], [0 < z < 1]. To 
observe the log dependence, we must switch to a system of ordinary differential equations 
after a certain value of f, beyond which the z-dependence has disappeared. 

The coefficient q is nevertheless obtained with good accuracy, since it is associated with 
a conservation law, obtained from the equation for the metric function A Eq.(p.3|), equivalent 
to the Smarr relation. The coefficient c^, in contrast, is error prone, since the asymptotic 
fall-off of the function B is numerically not well determined. 
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Figure 12. The mass M (a) and relative tension n (b) of the D = 5 and D = 6 nonuniform 

string branches arc shown versus the nonuniformity parameter A. The quantities are extracted from 
the first law and the Smarr relation, as well as from the asymptotic coefficients at f = 1. 

Reading off the values of the coefficient at f = 1 therefore leads to considerable dis- 
agreement of the values of the mass and tension with those values obtained from the first law 
and the Smarr relation. This is illustrated in Figure 12. 

Our final remarks concern the numerical method and the quality of the solutions. 

For the construction of the numerical solutions we use Newton-Raphson iteration. In 
each iteration step a correction to the initial guess configuration is computed. The maximum 
of the relative defect decreases by a factor of 20 from one iteration step to another. However, 
for large values of A convergence is slower. In this case we re-iterate the solution until the 
defect is small enough (about 10~^). Note, that this defect concerns the discretised equations. 
The estimates of the relative error of the solution (truncation error) are computed separately. 
They are of the order 0.001% for small A, but increase up to 1% close to the backbending 
point. On the second branch the errors first decrease with increasing A, but then increase 
again when the solutions become too steep at the waist. The errors also depend on the order 
of consistency of the method, i.e. on the order of the discretisation of derivatives. Usually 4**^ 
order gives reasonable results. For some solutions we used 6*^ order to check the consistency 
with the 4*^^ order solutions. 

We also monitored the violation of the constraints 

Ci = ^/fV^{G;-Gl), C2 = V^Gl. 

For small A the maximum of the constraints is less then 0.1, but it increases with increasing A. 
The maximum of the constraint C2 is smaller by a factor of 5 compared to Ci . On the second 
branch the maximum is of the order one. However, this large value appears in the vicinity 

of the waist (z = L/2 at the horizon), in a small region where the functions are extremely 
steep. Away from this region the violation of the contraints remains small. We also observed 



- 28 - 



that the condition dB/dr = at the horizon is violated at the waist when A becomes large. 
(A similar remark has been made in Q.) The violation of the condition dB/dr = at the 
waist might be related to the violation of the constraints. 

Increasing the number of grid points near z = L/2 yields smaller violation of the con- 
straints and the conditon dB/dr = at the horizon. Also, the violation of the constraint Ci 
is of the same magnitude as the maximum of the equations with the same weighting. 
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